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Capillary rise is one of the most well-known and vivid illustrations of capillarity, however, there 


is no complete solution has been obtained except some segmental solutions for asymptotic regimes. 


In this paper, we use singularity-free Bush’s equation and successfully obtain its series solution. 


The solution reveals that the capillary rise dynamics is mainly controlled by the Bond number 


and the Galileo number, of which the Bond number is a key parameter on the existence of the 


solution. Due to the poor rate of convergence of the series solution, we propose an approximate 


pga 


_( [82+ 
analytic solution h(t) = 229988 fı f H (v s) | at'| e ( 7E) i where a is a parameter. 


The proposed analytic solution has been verified numerically. 
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INTRODUCTION 


Capillary rise is one of the most well-known and vivid 
illustrations of capillarity. Knowledge of capillarity laws 
is important in oil recovery, civil engineering, dyeing of 
textile fabrics, ink printing, and a variety of other fields. 
It is the capillarity that brings water to the upper layer of 
soils, drives sap in plants, or lays the basis for operation 
of pens [1-25, 27-33, 35-38]. The capillary action has 
been studied by numbers of masters in science. Leonar- 
do da Vinci (1452-1519) recorded the effect in his notes 
and proposed that mountain streams may result from 
capillary rise through a fine network of cracks. Jacques 
Rohault (1620-1675): erroneously suggested that capil- 
lary rise is due to suppression of air circulation in nar- 
row tube and creation of a vacuum. Geovanni Borelli 
(1608-1675) demonstrated experimentally that h ~ 1/a. 
Geminiano Montanari (1633-87) attributed circulation in 
plants to capillary rise. Francis Hauksbee (1700s) con- 
ducted an extensive series of capillary rise experiments re- 
ported by Newton in his Opticks but was left unattribut- 
ed. James Jurin (1684-1750) independently confirmed 
h ~ 1/a; hence Jurin’s Law. Albert Einstein, whose first 
paper was submitted to Annalen der Physik in 1900, was 
on capillarity [1]. 

Washburn [3] developed the equation to describe the 
rate of penetration of liquids into small cylindrical cap- 
illaries based the laws of hydraulics. He assumed that 
the distance traveling in the initial turbulence period is 
negligible, and the Poiseuille region covers practically the 


entire flow. According to Poiseuille’s law, neglecting the 
air resistance, the rate of flow through a cylindrical takes 
the following form: Q = TAE (at + Lsa?), where a is the 
radius of the tube, h is the capillary penetration height, 


L, is the slip length. Q is the rate of flow and for a 
cylindrical tube it is given by Q = ra?44. The total 
effective pressure AP, consists of the external pressure 
Pe, the hydrostatic pressure P, and the capillary pres- 
sure given by Young-Laplace equation AP = 222288 
Combining the above relations and neglecting the s- 
lip length gives the so-called Lucas-Washburn equation: 
dh (Pa + Pp + 22298°). Tf there is no external pres- 
sure, a commonly used form of Lucas-Washburn equation 
is given by 8uh 44 +pgha?siny = 2accos6, where y is the 
incline angle of the tube with respect to the horizontal 


surface. 
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FIG. 1: The incline capillary rise 


When the gravitational force is negligible, a well- 
known Washburn’s law can be obtained with initial con- 
dition h(0) = 0: h = \/*%52°*¢, which predicts burst-like 
behavior with the velocity being infinite at zero time. 


This singularity of the solution highlights a deep inconsis- 
tency of the above equation [20]. Rideal [4] also obtained 
this equation through a rather simple derivation. 


Taking into account of the momentum of the liquid in 
the tube and end-effect drag on the fluid entering the 
tube, Brittin [6] derived a more rigorous formulations of 
mee Washburn equation as follows: h4% dh + 5(%)? + 

zaz uh 
viscosity, ø is the liquid- air surface tension 0 = yzyv, 0 is 


+gh= Zo cos, where p is density, u is the 


the wetting angle of the liquid, h is the height of capillary 
rise, a is the capillary radius, and g is the acceleration of 
gravity. This equation assumes the Poiseuille flow profile 
throughout the capillary. 


The most papular equation of capillary rise dynamics 


is hae + (By + ooo uh Fe +gh= Z0 608 0. Although 


this equation has been widely used in on literature, how- 


ever it has a singularity at t = 0, namely h = un — oo 
as t > 0, which would lead to an ill-posed problem. This 
equation can not even deal with the natural initial con- 
ditions h(0) = h(0) 
h(0) = \/20 cos 6/(pa) [5], neglecting such a logical draw- 


= 0. A formal remedy is to take 


back as the acceleration of liquid front at zero time being 
infinite, more criticism on this equation can be found in 
(8, 20, 31]. 

This singularity problem was firstly pointed out by 
Szekely et al.[8], who successfully removed the singulari- 
ty problem by composing the correct energy balance for 

2 
the entry flow: (h+ fa) GP +1.225(%)?4 Si Eh% +gh = 
r cos 6. Xiao et al.[28] considered the entrance pressure 


loss effects by using Dreyer’s model [13] and dynamic 
contact angle effect by Newman’s model [7] to further 
modified a Washburn equation: n +1. 028a) $} + 
0.958( 42)? + 4 (h +0.25a) È + gh = 22 cos 6. Bush [36] 
derived a similar equation with a diferenti in the coeff- 


cient of (4)?. 


ae 


For centuries of investigation on the capillary rise dy- 
namics, it is incredible that no complete solution has been 
obtained yet. Ones have to find some information for 
asymptotic regimes. How to match the different solu- 
tions between the regimes has never been solved and still 
remains as an unsolved problem. It would be natural 
thinking that ones can not satisfy those asymptotic solu- 
tions, and must be looking for a complete solution that 
can cover whole domain of time. In the following, we will 
mainly focus on Bush’s model and give its series solution 


and approximate analytic solution. 


To make the paper self-contained, the paper is organ- 
ised as follows. After this introduction, in Section 2 the 
Bush’s equation has been formulated. In Section 3 the 
series solution is obtained. In Section 4 an approximate 
analytic solution is proposed by using the initial condi- 
tions and the series solution. In Section 5 the numerical 
verification and discussions have been presented. Finally, 
section 6 concludes the paper. 


BUSH’S FORMULATIONS OF CAPILLARY RISE 
DYNAMICS 


In physics, the capillary dynamics process is a strug- 
gle between surface tension with the combining effect of 
wall resistance and gravity. In the initial stage, the sur- 
face tension plays a leading role, and later the combining 
effect of wall resistance and gravity takes over. 

Capillary rise dynamics is governing by a combina- 
tion of gravity, viscosity, fluid inertia and dynamics pres- 
sure. From continuum mechanics, the dynamics pro- 
cess must satisfy the conservation law of momentum: 
4 m(t) Z8] = F + f pv -v ndS, where the second 
term on the right-hand side is the total momentum flux 
Ta?ph? = mh, so the force balance on the column capil- 
lary may be expressed as 


(pra’h+ ma )h= 2naocos§ — pra?hg 


inertia added mass capillary force weight 
2 T 2 
= Ta ger — 2rah: T, 
—— 
—— 


g viscous force 
dynamical pressure 


(1) 


We can estimate the added mass m, from the change in 
kinetic energy as column rise from h to h + Ah, namely 
AE; = 4A(mU?), where m = me+mo+Meo (mass in the 
column, in the spherical cap and all other mass, respec- 
tively). In the column, me = ma?hp, u = U; in the spher- 
ical cap, Mmo = 2z a3 p; and in the outer region, radical 
flow extends to oo, nut u(r) decays. Volume preservation 
U/2. Con- 


tinuity gives 27a?u,(a) = 2rr?°ur(r) and leads to u,-(r) = 


requires: ta?U = 2ra?ur (a) leads to u,(a) = 


az u,(a) = £U . Thus the kinetic | in the far field 
is MU? = $f ur(r)?d $f. ulr 
3 


which leads to mao = Vora. . Hence the change in ki- 


)2 p2rrdr, 


201709.00105v2 


chinaXiv 


netic energy will be 


1 1 
OE, = zA Me + mo + Moo)U? + gm2U AU 
1 


1 
= zAV? + 5 (Me + mo + Meo)2UAU 


1 
= 5 Ala pAh)U* + (ma*hp + T a8 p)UAU. (2) 


7 
6 
For a typical capillary rise, a = 1074 ~ 5 x 1074m, u = 


The added mass is estimated as Ma = Mot+mM = trap. 
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FIG. 2: The dynamics of capillary rise 


1.5 x 10-°%8, p = 750kgm™* and top velocity can be 
at V = 0.04 ~ 0.778m/s, hence the typical capillary 
rise Reynolds number Re = eve = 10 ~ 195. At such 
low Reynolds number, if the tube is smooth one the flow 
in the tube is fully developed Poiseuille (laminar) flow, 
hence the viscous force is T, = —*#h, 
Combining all above relations, we have Bush equation 
of capillary rise dynamics [36]: 
7 dh 1 dh 


(h+ 69) Ga + ole 


8u , dh 
pa? dt 


4 +gh = a (3) 
pa 

and initial height and velocity boundary conditions 

h(0) = 0, h(0) = 0. 

Although no complete solution of Eq.(3) has never 
been obtained, some segmental solutions (as shown in 
Figure 3) for asymptotic regimes have been proposed by 
reduced form of Eq.(3). (1) Initial regime: h ~ 0, h ~ 0, 
infers A(t) = Peme hence h(t) = $ acot, (2) Quasi- 
steady stage (Lucas-Washburn equation), where the cap- 
illary force is compensated by gravity and viscous drag, 
the asymptotic solution corresponding to the short-time 
limit (t + 0) is given by A(t) = ,/ ®t. (3) Low- 
viscosity limit (Quéré equation), the capillary rise is giv- 
en by h = a (4) The capillary liquid will rise 
to stationary level, established by the balance of gravity 
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and capillarity, it gives H = ae 
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FIG. 3: Asymptotic regimes of capillary rise dynamics [36] 


However, those segmental solutions have not been 
matched into a single solution that valid for whole do- 


main of time. 


SOLUTION OF BUSH’S EQUATION OF 
CAPILLARY RISE DYNAMICS 


Eq. (3) is a nonlinear ordinary differential equation, up 
to now there is no closed form solution has been found- 
ed. Fortunately, it has an exact Taylor’s series solution 
h(t) = Xpo Cnt” as follows 


6o cos > 30 cos0(Tgpa? + 240 cos 0) 4 
E t t 
Tpa? 343p7a° 


864107 cos? 0 b 
= Te O(t 4 
1715p8a" oe) (4) 


The correctness of this solution has been verified by using 
symbolic mathematics software Maple. To understand 
the effect of surface tension, gravity and viscosity, let’s 
transform the Bush’s equation (3 into a dimensionless 


cos 0 dz 
8 1=0 
(\ Bo. =) ar 

(5) 
and the initial height and velocity boundary conditions 
are changed to 2(0) = 0, 4(0) = 0. In which, z = 4,7 = 
t — 2a0cosd@ = 
r H= pga ? T= 
represents the ratio of gravity and surface tension, and 


form as follows: 


dz i 7 Bo d@z i Te 
"Ir ' 12 cos 0 dr? ` 2`dr 


2 
= Bond number Bo = 22 


2 3 
Galileo number Ga = £? ma represents the ratio of gravity 
and viscosity. 
To have a fair quantitative understanding, we can 
rewrite the series solution in terms of Bo and Ga as fol- 


lows 
2(T) = 55 (40088) 2 ET s? 0)r4 
-= = GE cos? 0)r4 
~ cos? ( a) 7° 4+0(Bo~) (6) 


Qualitatively speaking, the Eq.(6) reveals that the Bond 
number Bo is playing a vital role in capillary rise dynam- 
ics, because it stands for the ratio of gravity and surface 
tension. In the very beginning (within the diffusion time 
a surface tension is mainly compete with gravity, on- 
ly after the diffusion time the viscous drag will come to 
against the surface tension. In other words, from the be- 
ginning of capillary rise, the Bond number Bo enters the 
process and later on the Galileo number Ga joins. 

The solution in Eq.(4) gives a important result, name- 


12 ø cos 
7 pa? ’ 


which later will be used in the proposal of approximate 
With the above solution in Eq.(4), 
there is no longer needs to divide the capillary dynamics 


ly, the non-zero initial acceleration A(0) = 
analytic solution. 


process into several asymptotic regimes, since the solu- 
tion in Eq.(4) is valid for the whore range of time do- 
main. From Eqs.(6,4) we can see that the Bond number 
Bo is really vital for the existence of the solution. For 
those equations without the added mass term fa mean- 
s that they have zero Bond number, which will lead to 
the infinite of height and velocity ! This is the source of 
singularity. 

The academic value of the series solution is enable us 
to see the importance of the added mass that can void 
the singularity. However, the rate of convergence of this 
solution is very slow and useless for practical applica- 
tion. It would be a natural attempts to propose a simple 


solution for the capillary rise dynamics. 


APPROXIMATE ANALYTICAL SOLUTION 


General speaking, the whole dynamics process can be 
qualitatively described as follows: in an infinite reser- 
voir, the capillary rise and velocity are zero in the initial 
state. Due to the effect of the surface tension, the cap- 
illary liquid obtains the initial acceleration (the initial 
acceleration must never be zero), begins to rise at a rel- 
atively uniform velocity, and the surface tension plays a 


dominant role in the ascending phase; however, as the 


capillary rises wall frictions and gravity begin to work 
in an attempt to prevent the rise of the capillaries, and 
their joint action succeeds in decelerating the capillaries 
at a point in time until the capillaries are finally stopped. 
Surface tension and wall resistance, gravity to achieve u- 
nity of opposites, capillary dynamics process are over, all 
attributed to calm. 

From the physical understanding of capillary rising 
phenomena, the capillary liquid will rise to stationary 
level, established by the balance of gravity and capillari- 
ty, the height and velocity will tends toh — H = Ae cost 
and h + 0 at t > ov, respectively. To satisfy the con- 
ditions, we can propose that the capillary rise takes the 


following form 
h(t) = H[1 — fe), (7) 


where the exponent 8 and function f(t) are to be deter- 
mined. 

Substitute the Eq.(7) and h(t) = —Hf’(t)e7® + 
H£f(t)e~® into the initial conditions h(0) = 0, h(0) = 
0, the initial condition for the function f(t) is given by 


f(0) = 1, f'(0) = £. (8) 


There are many possibility on the construction of f(t) 
which can satisfy the condition in Eq.(8), for instance, 
f(t) = 1 + 6t+ at’, in which the exponent d > 2 and 
a is a constant. Since capillary rise dynamics process 
usually takes a millisecond, so the lower order of time t 
will dominate the process, based on this understanding, 
we set the function f(t) = 1+ St + atè. Hence, the 


capillary rise can be proposed as follows 
A(t) = H[1 — (1+ Bt + at?)e"™*]. (9) 


Thus, we have capillary velocity h(t) = H(aft? — 3at + 
B?)te-** and acceleration h(t) = —H(a6?t® + Bt — 
6aft? + bat — B?)e—**. 

The capillary rise in Eq.(9) satisfies both initial con- 
dition and infinity condition. The unknown constants 8 
can be determined by the exact series solution in Eq.(4), 
t 1p coos? H . Hence, 
by equaling h(0) = H8?, we obtain the exponent 


12 ocos 6 
B= Le 4/2. (10) 
7H pga Ta 


To find the parameter a, we can expand the capillary 


which gives the acceleration h(0) = 


rise in Eq.(9) in Taylor series, and get its coefficient of 


the fifth order team of t° as — 3H6?(8? + 15a), and let 
it equal to the coefficient of të in Eq.(4), namely 


— 5 H167(B° + 18a) = ae et (11) 
which gives 
Tuocos@ 1 ( iN) 
sat (yE) aa 


It is clear to see that œ is crucial parameter which links 
to all constants the capillary rise dynamics. It is worth 
to mention an interesting thing that a can be set to zero 
by adjusting the radius. 
Finally, we obtain the capillary rise h(t) 
_ 20 cos 


6g 3 
h(t —— ;1-]i —|t t 
= 6 pga i ( 2) t i 


and the capillary velocity h(t) 


| 6 _(,/88 
(e 8 at + s) te (vat, 


(14) 


We) 
(13) 


dh(t) _ 20 cos0 


dt pga 


as well as the capillary acceleration h(t). 

The above solutions in Eqs.(13) and (14) are valid for 
the whole domain of time t € [0, 00), for my knowledge, 
which have never been proposed before. 

The solution in Eq.(13) reveals that the capillary rise 
h(t) is mainly controlled by the H = accost with decay 


_(,/ 8a 
rate of e (Vaat where the radius a is the only dominate 
parameter, the smaller radius the fast decay. 


© NUMERICAL VERIFICATION AND 
DISCUSSIONS 


Verification of the solution in Eq.(13) 


To verify the correctness of the solution in Eq.(13), we 
can compare the result obtained numerically by Zhmud 
et al.[20]. for comparison study, the data of the diethyl 
either in glass shown in Table I. 


TABLE I: The diethyl ether in glass with arbitrary radius 


[EZ] o [E9] 0 g|5] e(=S] 
2.2.1074 1.67-10-2 26° 9.81 710 


In this case, for any radius a, we have the capillary rise 
H= 43L - 1076, our analytical capillary rise is given by 


4.31 - 1076 
h(t) = ———_{1- 
2.9 3.821- 1071? 1.626 _ 2.9 
[+ at bf a “aT yt Je" va}. (15) 


If we set the radius a = 0.0005, we have the capillary rise 
as follows 


A(t) = 0.00862 [1 — (1 + 129.68¢ — 23112.54t%)e7 179") , 

(16) 
which is plotted in the Figure 4. From this figure we can 
see that our solution agrees with [20] very well. 
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FIG. 4: Comparison study: a = 0.0005m and @ = 26° 


The capillary rise velocity is given by 


h(t) = —25836¢(t + 0.06)(t — 0.08)e~ 129", (17) 


which is plotted in the Figure 5. The capillary rise accel- 


eration is given by 


h(t) = 3350496(t + 0.006) (t — 0.008) (t — 0.096)e~ 179-68", 
(18) 
which is plotted in the Figure 6. 


Influence of radius change 


If keeping the wetting angle as 0 = 26° [20], and just 
adjusting the radius of the glass in Eq.(15). We have 
graphs in the Figure 7, which reveals that capillary rise 
is very sensitive to the radius. The smaller the radius 
in which the liquid can travel, the further up it goes. 


0.55 


0.444 


—m— present 


capillary rise velocity 


FIG. 5: Velocity profile for a = 0.0005m and 6 = 26° 
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FIG. 6: Acceleration profile for a = 0.0005m and 6 = 26° 


Therefore, the glass tube radius is a crucial parameter 
for capillary rise dynamics. The radius has also strong 
influence on the capillary velocity in the Figure 8 and 
acceleration which can been seen in the Figure 9. The 
smaller the radius the fast the capillary goes. 


Discussions 


If we put all the pictures together, we can have a fair 
picture of the capillary rise dynamics as follows: the cap- 
illary rise onsets from zero initial height and velocity, at 
t = 0 as the first driving force, the surface tension sup- 
plies a kick-off acceleration in the amount of H?, this 
surface tension steadily drives the capillary rising at an 


almost uniform speed until reaches its peak at certain 


capillary rise 


capillary rise velocity 


FIG. 8: The capillary rise velocity for four different radius 


point, then the capillary speed is gradually reduced to 
zero due to the combining effect of wall viscous friction 
of galls and gravity. In the same time, the capillary accel- 
eration decrease and reaches its lowest point before the 
increase to its final static state. 


CONCLUSIONS 


Since the capillary rise generally has low Reynolds 
number, the flow in the tube can be considered as a lam- 
inar flow. To avoid the singularity problem, the added 
mass must appear in the governing equation. From the 
dimensionless form of Bush’s equation, we found that the 
capillary rise dynamics is mainly controlled by the Bond 
number and the Galileo number. We have successfully 


capillary rise acceleration 


FIG. 9: The capillary rise acceleration for four different radius 


obtained the perturbation and series solutions. Due to 
the poor rate of convergence of this solution, we pro- 
posed a simple approximate analytic solution by using 


series solution and initial conditions. We have verified 


the proposed analytic solution by some numerical exam- 


ples. 


Last but at least, it might be worth to point out that 


all previous solutions are for the non-oscillatory regime, 


however, when the liquid surface reaches h = H, the sur- 
face will oscillating for while before asymptotically stop. 
Within the oscillatory regime, the above solutions should 
be modified a little bit, which will be investigated by 
Homotopy Analysis Method [26, 34] in the forthcoming 
dg paper. 
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